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Abstract. We propose a simple graphical approach to steady state solutions of the cable equation 
for a general model of dendritic tree with tapering. A simple case of transient solutions is also briefly 
discussed. 



1. Introduction 
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The function of many physiological systems depends on branched structures that exist both at 
the tissue (e.g. nervous plexi, lungs, and the vascular and lymphatic systems) and the cellular level 
(e.g. neurons). Of particular interest, local and global propagation of electrical signals within the 
Q -1 ! nervous system depends on the integration, processing, and further generation of electrical pulses 
that travel through neurons. In turn, the tree-like morphology of neurons facilitates simultaneous 

t-h ■ signaling to cells located in different places and over long distances. 

> ; 

Neuronal morphology is typically modeled by assuming that the shape of small neuronal segments, 
£^ \ or neurites, is approximated by cylinders of different diameters. As a consequence, cable theory [25] , 
[46], [UJ, [IB], [IS] j [50] plays a central role in the theoretical and experimental study of electrical 
conduction in neurons; see, for example |], [II], [5H], [H], [12], [2U], [21], [22], [23], [IS], [S3], [57J 
and references therein. Notably, one of the most interesting results from recent theoretical work is 
that geometrical properties of neuronal membranes may exert powerful effects on signal propagation 
■ even in the presence of voltage-dependent channels [60] . 



^ . Theoretical research involving realistic neuronal morphologies is typically done by numerically 
solving systems of cable equations defined on cylinders with different radii, and assuming that 
voltage and current are continuous functions of space and time. To the best of our knowledge, 
graphical methods seem have not been widely applied yet in the mathematical modeling of neurons. 
Graphical methods are very useful and popular in different branches of modern physics. It is worth 
noting, for example, Feynman diagrams in quantum mechanical or statistical field theory [3], [I], 
[TO] . [25] . [26] . [27] . [32] , [37] . [63] . Vilenkin-Kuznetsov-Smorodinskii approach to solutions of n- 
dimensional Laplace equation [38], [12], [51], [55], applications in solid-state theory, etc. A goal of 
this paper is to make a modest step in this direction (see also [JJ, [U] and references therein). We use 
explicit solutions from recent papers on variable quadratic Hamiltonians in nonrelativistic quantum 
mechanics [12], [15], [TB], [17], [JS], [35], [3S], [55], [5H] to describe steady state and transient 
solutions to linear cable equations modeling neurites with non-necessarily constant radius. 
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2. Cable Equation with Varying Radius 

At a closer view, neurites can be regarded as volumes of revolution, defined by rotating a smooth 
function r = r(x) representing the local radius of the neurite where x represents distance along the 
neurite. As a result, the cable theory implies the following set of equations [3l] . [4"T] : 

2 I ^ S ^ ^ S 1 + ( ( ^ r \ 2 (2 1) 

m dx dx 1 dx y \dx ) 

V dV 

Im= R~ + Cm ^t> (2 ' 2) 

Here, represents the voltage difference across the membrane (interior minus exterior) as a devi- 
ation from its resting value, I m is the membrane current density, I = I a is the total axial current, 
R m is the membrane resistance, Ri is the intercellular resistivity and C m is membrane capacitance 
(more details can be found in [31], [32] and [UJ). Differentiating equation (12. 3p with respect to x 
and substituting the result into (12. ip with the help of ( 12. 2 p one gets 

dV\ i , (dr\ 2 R m 1 d ( 2 dV 



which is the cable equation with tapering for a single branch of dendritic tree (see [H], [31], [17] 
and [57] for more details). 

We shall be particularly interested in solutions of the cable equation (12 .4p corresponding to 
termination with a "sealed end", namely, when at the end point x — x% the membrane cylinder 
is sealed with a disk composed of the same membrane. In this case, the corresponding boundary 
condition can be derived by setting 

h = Kr 2 I m , (2.5) 
at x — x\. Then, in view of ( I2.2p -( F2T3|) . one gets [45J: 

= 0, t > 0. (2.6) 



= 0, t > 0, (2.7) 

X=XQ 



v + r r dv + Rmdv 

In a similar fashion, at the somatic end one gets 

where R s is the somatic resistance and C s is the somatic capacitance [20]. We shall use these 
conditions for the steady-state and transient solutions of the cable equation. (Later we may impose 
similar boundary conditions at the points of branching.) 

In this Letter, we shall first concentrate on steady-state solutions of the cable equation, when 
dV/dt = 0. Then 

R nl (1 f r 2dV\ (xo < a; < ;Ei) (2 . 8) 




IrRi dx V dx 
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and 

^Un = V 0t (^z + B (x) V 



dx 



= 0, B( Xl ) = ^. (2.9) 



This boundary value problem can be conveniently solved (by a direct substitution for each branch of 
the dendritic tree) in terms of standard solutions of this second order ordinary differential equation 
as follows 

\r ( \ w( v C (x) + B (xi) S (x) 

"'^"WcW + JWsW (2 ' 10) 

where C (x) and S (x) are two linearly independent solutions of the stationary cable equation (12.81) 
that satisfy special boundary conditions C (xi) = 1, C (xi) = and S (x±) = 0, S' (xi) = —1. Then 

dV 

— + B(x)V = (x <x<xi) (2.11) 
ax 

with the corresponding current density/voltage ratio function B (x) given by 

B(x) V ' C (x) + B ( Xl ) S' (x) 

[) ~ V~ C(x) + B( Xl )S(x) [ ' 

in term of the standard solutions C (x) and S (x) . Throughout this Letter, we shall refer to a case, 
when B (x) > (x < x < x\) , as the case of weak tapering. An opposite situation, when B (£) = 
at certain point Xq < £ < x\ and an inverse of the current may occur, shall be called a case of the 
strong tapering. (A case of strong tapering has been numerically discovered in [28].) 



3. Tapering with Analytic, Asymptotic and/or Numerical Solutions 

In this Letter, we consider a general model of a dendrite as a (binary) directed tree (from the 
soma to its terminal ends) consisting of axially symmetric branches with the following types of 
tapering. 



3.1. Cylinder. Here, r = r =constant, < x < L and the cable equation takes the simplest form 

d 2 V A 2 = 

dx 2 ' 2R { 



X 2 Pj = V, A 2 = ^ (3.1) 



with a familiar solution [31], [15], [%T] : 



V( \ \r cosh (( L ~ x ) / A ) + Ajg ^ sinh (( L ~ x ) / A ) /oo^ 

{ 1 cosh (L/A) + \B L sinh (L/A) 1 J 



subject to boundary conditions 



V(0) = V , B L V(L) + ^(L) = 0. (3.3) 

dx 



Then 

A J B L + tanh((L-a;) /A) 
A + A 2 S L tanh((L-x) /A) 
(See [31], [17], [19], [50], [59] and references therein for more details.) 



(3-4) 
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3.2. Frustum (Cone). Here, r = r (x) = r + Tl L r ° x with < x < L. The steady-state solution of 
the corresponding cable equation 

A 1/4 



r ^ +2 ^ = T^' A4= v^: — — (3 - 5) 



L 



subject to boundary conditions (13. 3 j) are given by [7] 



Here, the standard solutions that satisfy S (0) = 0, S" (0) = 1 and C (0) = 1, C" (0) = can be 
constructed as follows 

2Lr 3/2 

S (x) = - ° ^_ (ji^) h (tiy/Pj - Ji (//VFb) #i (/V^] (3.7) 

and 

C(x) = [/x v ^iro(^) + 2ir 1 (^)] A /^/ 1 (/x v / ^) (3.8) 

V ri 

in terms of modified Bessel functions I v (z) and K v (z) of orders v = 0, 1 (different aspects of the 
advanced theory of Bessel functions can be found in [2J, [5], [6], [21], [H], [13], [H], [61] and [62]). 

3.3. Hyperbola. If r = a cosh on an interval xo < x < xi, the cable equation (12. 4 ft takes the 
form 



1/ + C m R m —— cosh 



at / V a / 2i? i 



a; — fe\ 91/ fx — b\ d 2 V 

2 sinh — — h a cosh 



dx \ a J dx 



(3.9) 



This special case of tapering is integrable in terms of elementary functions [30] (see also [12] and [17] 

for a similar problem related to a model of the dumped quantum oscillator). For the steady-state 
solutions one obtains the following equation 

V" + 2Xtanh(Xx + 5)V = i4V (3.10) 

with new parameters 

A = i S = -Xb, l4 = ^w~- (3.11) 
The corresponding two linearly independent solutions, namely, 

T , / s sinh (ux + 7) . . cosh (ux + 7) ~n" 

^W = oosh(Ax + lr ^W= co S h(A, + I) - " + (« 2 > 

can be verified by a direct substitution for an arbitrary parameter 7. 
The required steady-state solution of the boundary value problem 

V(x ) = V Q , BV(x 1 ) + ^(x 1 ) = (3.13) 



is given by 



where 



and 
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VM-y C(x 1 -x) + BS{x 1 -x) 

v {x) ~ Vo c(x 1 -x ) + BS(x 1 -x y (3 ' 14) 

C( Xl - x) ^o sH X( b -^) C -^^§ (3.15) 
. i f \ / 1 w Asinh(/i(xi — x)) 



(See [30] for more details.) 

3.4. General Case of Axial Symmetry. One can use numerical methods and/or WKB-type 
approximation in order to obtain standard solutions. For example [4"T] . 

V « 1 [Ae^ x) + Be~^ x) ] , (3.17) 

^Jr 2 {x)p{x) 



where 



2 ds Ri x 1/2 



P{X)= ( r^iT' ' ' P(t) (3 - 18) 



(See [31] and [17] for further details.) 

4. A Graphical Approach 

Graphical rules for steady-state voltages and currents in a model of dendritic tree with tapering 
are as follows. 

4.1. Single Axially Symmetric Branch with Arbitrary Tapering. For a single branch with 
tapering voltage and current density/voltage ratio are given by 

V(x) = C(x 1 -x) + B(x 1 )S(x 1 -x) 

[ ' C ( Xl -x ) + B ( Xl ) S ( Xl - x ) [ 0) [ } 

= (C (xi — x) + B (xi) S (xi — x)) V (xi) , 
p ,s _ C (xi -x) + B (xQ g (xi - x) _ V (x) 
lJ C(x 1 -x) + J B(x 1 )5(x 1 -x) V(x)' l ' J 
respectively (see Figure 1). 



Figure 1. A Single Branch with Tapering. 
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4.2. Junction of Three Branches with Different Types of Tapering. The internal potential 
and current are assumed to be continuous at all dendritic branch points and at the soma-dendritic 
junction [45]. We consider a general case when each branch has its own tapering, say r = r (x) , 
T\ = ri (x) and r 2 = r 2 (x) (see Figure 2). Then 

V (x 12 ) = V 1 (x 12 ) = V 2 (x 12 ) , (4.3) 

r 2 (x 12 ) B (x 12 ) = r\ (x 12 ) B x (x 12 ) + r 2 (x 12 ) B 2 (x 12 ) . (4.4) 




Figure 2. A Junction of Three Different Branches. 
The total ratio constant B (xi 2 ) at the branching point x± 2 is given by the following expression 



B(x 12 ) = r 4^B 1 (x 12 ) + r 4^\B 2 (x 12 ) 
r 2 (x 12 ) r l (x 12 ) 

= B(B 1 (x 1 ),B 2 (x 2 )) 

rj (x 12 ) C[ (xi - x 12 ) + B 1 (xi) S[ (xi - x 12 ) 

r 2 (x 12 ) Ci {xx - x 12 ) + B x (xi) S x [x x - x 12 ) 

, f 2 2 (^12) C' 2 (x 2 - x 12 ) + B 2 (x 2 ) S' 2 (x 2 - x 12 ) 



r 2 (x 12 ) C 2 (x 2 - x 12 ) + B 2 (x 2 ) S 2 (x 2 - x 12 ) 



Then the ratio constant B (xq) is 

B (X ) : 



C (X 12 - Xq) + B (x 12 ) S' (X 12 - Xq) 
C (X 12 - Xq) + B (x 12 ) S {X X2 - Xq) 

with the coefficient B (xi 2 ) found by the previous formula (14. 5p . 



(4.5) 



(4.6) 



4.3. Junction of (n + 1)-Branches. In a similar fashion, at the branching point x a , one gets 

V (x a ) = V, (x a ) = V 2 (x a ) = ... = V n (x a ) , (4.7) 



r 2 (x a ) B (x a ) = r 2 1 (x a )B 1 (x a ) +rl(x a )B 2 (x a ) + ... + r 2 (x a )B n (x a ) 

n 

= r] (x a ) Bi (x a ) . 



(4.8) 
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(See Figure 3.) 



V 




Figure 3. A Junction of (n + 1)-Branches. 



Combination of the above graphical rules results in a simple algorithm of evaluation of voltages 
and currents in the model of dendritic tree under consideration as follows. 

Evaluate constants B (x a ) for all branching points of the tree: (a) first apply formula (I4.5P for all 
open notes; (b) remove the above nodes from the tree and keep repeating the previous step until 
you reach the root of tree (soma). 

In order to find voltage at a point x of the dendritic tree, follow the path xq — > x and multiply 
the initial voltage V (xq) by consecutive corresponding factors from formula (14.11) changing at each 
intersection of the tree. The ratio of voltages V (x a ) and V {xp) at two terminal points, can be 
determine in a graphic form by the previous rule applied to the shortest path x@ — > x@. 



Our formulas fl4.5p - fl4.6p define ratio coefficients B (x a ) for all vertexes for the standard node on 
Figure 2. For the corresponding voltages, one can write 



5. Examples 



V(x ) 



[C (xi2 -x ) + B (x 12 ) S (x 12 - x )) V (x 12 ) 

[C (x 12 - Xq) + B (x 12 ) S (X 12 - X )] V (X 12 ) 

x [C (xi - xi 2 ) + B (xi) S (xi - X12)} V (x\) 

[C (X l2 ~X )+ B (x 12 ) S {x 12 - Xq)} V (x 12 ) 

x [C (x 2 - xi 2 ) + B (x 2 ) S (x 2 - X12)] V (x 2 ) 



(5.1) 



and 



V (xi) _ C (x 2 - x 12 ) + B (x 2 ) S (x 2 - x 12 ) 

V (x 2 ) C (xi - X12) + B (xi) S (xi - xn) 



(5.2) 



Further examples are left to the reader. 
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6. Transient Solutions 



6.1. A Single Branch with Smooth Tapering. Let us consider the cable equation (12. 4p for 
a single branch with an arbitrary smooth tapering r = r (x) on the interval Xq < x < x\. The 
separation of variables 



V(x,t) = e -( 1+a2 > /Tm U (x) , r r 



in results in 



ld_ 

r dx 



,dU 
dx 



+ 0^17 = 0, 
dx 



to 



" - -a , 



Rr. 



(6.1) 



(6.2) 



where a is a separation constant. The boundary condition at the sealed end (12. 6ft takes the form 

= 0. (6.3) 



\dx 2 



X = X\ 



A general solution of this problem can be conveniently written (for each branch of the dendritic 
tree) as follows 

C (x, uj) + -uj 2 S (x, oj) 



U(x) = U(x,u) = A 



(6.4) 



where A is a constant and C (x, u) and S (x, u) are two linearly independent standard solutions of 
equation ( 16. 2 j) that satisfy special boundary conditions C (x%,cj) = 1, C (xi,u) = and S (xx,co) = 
0, S' (xi,uj) = 1. Then the boundary condition (12. 7ft at the somatic end x = Xo, namely, 



dU 

dx 



+ 



T 



R. \ r. 



C 



2Q 



s 2 



U 



o, 



X=XQ 



results in a transcendental equation 



1-^-11 

I'm 



R 



2R, 



m 2 
CO 



Ri 
R x 



C (x ,u) + \u 2 S' {x ,uj) 
C (x , uo) + \uj 2 S (x , oj) 



(6.5) 



(6.6) 



for the eigenvalues oj. (There are infinitely many discrete eigenvalues [13] and [29], Ince?.) The 
corresponding eigenf unctions U n = U (x, co n ) = A n u n (x) are orthogonal 



("mi fimn \Uni ^n) 

with respect to an inner product that is given in terms of the Lebesgue-Stieltjes integral 
also Appendix and [33], [51] and [52]): 



(u,v) : 



u (x) v (x) r (x) ds 



(6.7) 
(see 

(6.8) 



1 C 

+ -r 2 (xi) u (xi) v (xi) + 7^h r2 ( x o) u (x ) v (x ) . 



A formal solution of the corresponding initial value problem takes the form 
V(x,t) = V(x,oo) 



exp 



-, . Rm 2 

1 + W^ 



Tr. 



(6.9) 
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where V (x, oo) is the steady-state solution, u = u n are roots of the transcendental equation (16.61) 
and the corresponding eigenf unctions are given by 

1 



U n (x) = C (X, U n ) + -U n S (x, U n ) . 



(6.10) 

Coefficients A n can be obtained by methods of Refs. [13] and [20] with the help of the modified 
orthogonality relation (16. 7p as follows 

(V(x,0)-V(x,oo),u n (x)) 

[U n (x),U n [X)) 

Substitution of (16. lip into ( \6.9\\ and changing the order of summation and integration result in 

V(x,t) = V{x,oo) (6.12) 

G(x,y,t)(y(y,0)-V(y,oa)) dfjt(y), 



+ 



Supp a 



where 



G (x, y,t) = ^2 ex P 



u n (x) u n (y) 



Ur 



(6.13) 



is an analog of the heat kernel. Infinite speed of propagation. Method of images. 

6.2. A Single Branch with Piecewise Typering. In the case of a piecewise tapering 

r (x), x < x < x i 

T\ (x) , Xoi < X < X\ 

with r (^oi) = Ti (xqA , in a similar fashion, one can write 

A [C (x, u) - (f - l) + ^uA S (x, u] 
A x [Ci (x,u) + \u 2 Si (x,u)\, 

provided Cq{xq,uj) = 1, C (xo,u) = and So(xo,u) = 0, S' q (xq,u) = 1 and Cx(xx,u] 
C[ (xi,lu) = and Si (xi,u) = 0, S[ (xx,u) = 1. Continuity and smoothness of the solution 
point Xqi, namely, 

V (x 01 ) _ V (x+) 



U (x, u) 



, Xq < X < X 1 
Xoi < X < X\ 



(x q1 ,oj) + \<u 2 S[ (x x,u) 



U (x 01 ) U (*+) ' 
results in the following equation for the eigenvalues 

C> (x 01 ,u) - (f (g - l) + ^ ) g (wj c[ 

C (x 01 ,u) - (f - l) + 2c^ 2 ) fcoi, w) Cl ( Xo1 ' w ) + ^ 2Sl ( Xo1 ' w ) ' 
Introducing 

' - (0) (x)M 0) (x 01 ), * <*<* 01 



U n (X 



a 



(i) 



where 



«1 0) (x 



(x) = Cq (X, U n ) - 



{x) /Un (xqx) , X i<X<Xi 

Ri 



(6.14) 



(6.15) 

= 1, 
at the 

(6.16) 



(6.17) 



(6.18) 



s \ 1 m 



1 I + TTFT^n ) S (X, U n ) 
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U$P (x) = d (x, UJ n ) + -jU 2 Sx (X, U n ) , 

one can obtain a formal solution in the form fl6.12p - fl6.13p once again. Further details are left to 
the reader. 

7. Summary 

In this Letter, we propose a simple graphical approach to steady state solutions of the cable 
equation for a general model of dendritic tree with tapering. A simple case of transient solutions is 
also briefly discussed. 
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Appendix A. Modified Orthogonality Relation 



We consider the Sturm-Liouville type problem, 

Lu + Xpu = 0, 

for the second order differential operator 



Lu = — 

ax 



ax 



q fx) u, 



(A.l) 



(A.2) 



where k, q and p are continuous real- valued functions on an interval [xo,Xi] , k and p are positive 
in [x ,Xi] , k! exists and is continuous in [x ,Xi] , subject to modified boundary conditions 

v! (x ) + (a + b X) u (x ) = 0, (A. 3) 

v! (xx) + (ai — b±\) u (xi) = 0, 

where ao,bo > and a\, b\ > are constants. With the help of the second Green's formula (see, for 
example, [3D]), 

r\ v Lu-uL v) i x=k L^-u d A X \ (A.4) 
Jxo V dx dxj Xo 

for two eigenfunctions u and v corresponding to different eigenvalues 

Lu + Xpu = 0, Lv + ppv = 0, A 7^ p (A. 5) 

one gets the following orthogonality relation [13J: 

(A.6) 



u (x) v (x) pdx + b\k (xi) u (xi) v (xi) + bok (xo) u (xq) v (xq) = 0. 



^■0 



Here, the modified inner product 



(u,v) 



uv dp 



(A.7) 



Supp fj, 
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-x 1 

u (x) v (x) pdx + b\k (xi) u (xx) v (xi) + bok (xq) u (xq) v (xq) 



■CO 



is defined in terms of the Lebesgue-Stieltjes integral [31]. The modified orthogonality relation (1A.6I) 
holds also in the case of a piecewice continuous derivative k' on the interval [xo, xi] . 

The junction of three branches (see Figure 2) can be considered in a similar fashion. Suppose 
that 



d 

LiUi + XpiUi = 0, LiU = — 



h {X) fx 



-qi{x)u (A.8) 



T 



with k — 0, 1, 2 for three corresponding branches, respectively, and boundary conditions are given 
by 

u' (x ) + (o + b X) u (x ) = 0, (A.9) 
u' (xi) + (ai — bxX) u (xx) = 0, 
u' (x 2 ) + (a 2 - b 2 X) u (x 2 ) = 
at the terminal ends. Introducing integration over the whole tree T by additivity, 

(vLu — uLv) dx = {vqLqUq — uqLqVq) dx (A. 10) 

T Jx 

PXl PX2 

+ / {vxLxUx - UiLtVx) dx + (v 2 L 2 u 2 - u 2 L 2 v 2 ) dx, 

J X\2 J X\2 

and applying the Green formula (1A.4I) for each branch, one gets 

(vLu - uLv) dx = k (x 12 ) (v (x 12 ) u' (x 12 ) - u (xi 2 ) v' (x 12 )) (A. 11) 

-ki (X 12 ) (v i {X 12 ) u'x (X 12 ) - Mi (x 12 ) v[ (x 12 )) 
-k 2 (x 12 ) (v 2 (x 12 ) u 2 (x 12 ) - u 2 (x 12 ) v' 2 (x 12 )) 
-k (x ) (v (xq) u' (x ) - u (x ) v' (x )) 
+ki (x{) (v i (xi) u[ (Xi) - Mi (xi) v[ (Xi)) 
+k 2 (x 2 ) {v 2 (x 2 ) u 2 (x 2 ) - u 2 (x 2 ) v' 2 (x 2 )) . 

We shall assume that the following continuity conditions: 

u Q (xi 2 ) = «i (xi 2 ) = u 2 {x X2 ) , (A. 12) 

h (x i2 ) u' (x i2 ) = ki (x i2 ) u\ (x i2 ) = k 2 (x i2 ) u' 2 (x 12 ) 

hold at the branching point xi 2 . In view of of the boundary conditions (IA.9I) . the modified orthog- 
onality relation takes the form 

X 1 Ml PX 1 

u (x) v (x) pdx + / u(x)v(x) pdx + / u(x)v{x) pdx (A. 13) 

XO J XQ Jxo 

+b k (x ) u (x ) v (x ) + bik (xi) u (xi) v (x x ) + b 2 k (x 2 ) u (x 2 ) v (x 2 ) = 0. 

The case of junction of (n + l)-branches (see Figure 3)is similar. In general, for an arbitrary tree, 
one may conclude that only the terminal ends shall add additional mass points to the measure, if 
the corresponding boundary and continuity conditions hold. Further details are left to the reader. 
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